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ABSTRACT 

We present the analysis of the time series observations of C0R0T 102699796 ob- 
+3 ■ tained by the CoRoT satellite that show the presence of five independent oscilla- 

tion frequencies in the range 3.6-5 c/d. Using spectra acquired with FLAMES@VLT, 
we derive the following stellar parameters: spectral type F1V, T e ff=7000±200 K, 
log(g)=3.8 ± 0.4, [M/H]=-l.l±0.2, usim=50 ± 5 km/ a, L/L =21+^. Thus, for the 
first time we report the existence of a metal poor, intermediate-mass PMS pulsating 
star. Ground-based and satellite data are used to derive the spectral energy distribu- 
tion of CoRoT 102699796 extending fr om the optical to mid-infrared wavelengths. The 
SED shows a significant IR excess at wavelengths greater than ~ 5ii. We conclude that 
CoRoT 102699796 is a young Herbig Ae (FIVe) star with a transitional disk, likely 
associated to the HII region [FT96]213.1-2.2. 

The pulsation frequencies have been interpreted in the light of the non-radial 
pulsation theory, using the LOSC code in conjunction with static and rotational evo- 
lutionary tracks. A minimization algorithm was used to find the best-fit model with 
M=1.84 M , T off =6900 K which imply an isochronal age of t~2.5 Myr. This result 
is based on the interpretation of the detected frequencies as g-modes of low-moderate 
n-value. To our knowledge, this is the first time that such modes are identified in a 
intermediate-mass PMS pulsating star. Since CoRoT 102699796 lies in the region of 
^ | the HR diagram where the S Set and 7 Dor instability strips intersect, we argue that 

the observed pulsation characteristics are intermediate between these classes of vari- 
ables, i.e. CoRoT 102699796 is likely the first PMS hybrid 7 Dor-6 Set pulsator ever 
studied. 

Key words: stars: pre-main-sequence - stars: variables: T Tauri, Herbig Ae/Be - 
stars: variables: S Scuti - stars: fundamental parameters - stars: abundances - infrared: 
stars. 
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* The CoRoT space mission, launched on December 27th 2006, 
has been developed and is operated by CNES, with the contribu- 
tion of Austria, Belgium, Brazil , ESA (RSSD and Science Pro- 
gramme), Germany and Spain. 

f Based on observations collected at the European Southern Ob- 
servatory, Paranal, Chile. Proposal n. 082.D-0839A 
| E-mail: ripepi@oacn.inaf.it 



1 INTRODUCTION 

Asteroseismology of Herbig Ae stars offers a unique means 
to probe their interior structure and to compare to evolu- 
tionary models. It is now well established that these stars 
during contraction towards the main sequence cross the 
instability strip of more evolved stars. These young, pul- 
sating intermediate-mass stars are collectively called PMS 
5 Set and their variability is characterized by short pe- 
riods (~30m-h5 h) and small amplitudes (from less than 
a millimag to a few hundredths of magnitude, see, e.g. 
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Kurtz fc Marangj 1 19951 ; ICatalal 120031 ; iRipepi et all 120061 ; 
Zwintj|2008t ). 



The first theoretical investigation of the PMS instability 
strip based on non linear convective hydrod ynamical models 
was carried out bv lMarconi fe Pallal |l998l) . who calculated 
its topology for the first three radial modes. These authors 
also found that the interior structure of PMS stars cross- 
ing the instability strip is significantly different from that of 
more evolved Main Sequence stars (with the same mass and 
temperature), even though the envelope structures a r e sim - 
ilar. The subsequent theoretical work bv lSuran et alJ l|200ll ) 
made a comparative study of the seismology of a 1.8 PMS 
and post-MS star. They found that the unstable frequency 
range is roughly the same for PMS and post-MS stars, but 
that some non-radial (g) modes are very sensitive to the deep 
internal structure. More recently, iGrigahcene et all (|2006l ) 
produced a theoretical insta bility strip for PMS stars for 
the first seven radial modes; IRuoppo et al.l (|2007l ) derived 
a model based relation between the large frequency sepa- 
ration and the stellar luminosity and effective temperature 
and developed a tool to com pare theory and observatio ns in 
the echelle diagram. Finally |Pi Criscienzo et al.l (|200fii ) ap- 
plied the ATON evolutionary code to the computation of 
detailed grids of standard (non-rotating) and rotating pre- 
main sequence (PMS) models and computed their adiabatic 
oscillation spectra. 

From th e observational point of view, multi-site cam- 
paigns (e.g., IRipepi et al.l 120031 ; iBernabei et al.l 120091') and 
space observations with MOST (e.g., Zwintz et al. 20091 . for 
NGC2264) provided us with good sets of frequencies to be 
compared wit h the models. In addit i on, spectroscopic s tud- 
ies (see, e.g. iFumel fc Bohml |200S| ; iBohm et al l 120091 for 
the cases of HD 104237 and RS Cha, respectively) demon- 
strated to be very valuable for line profile an alysis and direct 
mode determination usin g tools like F2D (Kcnncllvl 1 19941 ; 
iKennellv fc Walker! Il996l ) or FAMIAS (|Zimal l200Sft . How- 
ever, the full exploitation of the PMS 8 Set instability strip 
is still far from being realized. Dedicated observations of 
a large sample of PMS stars covering the whole instability 
strip are needed to accomplish this goal. 

In this context, a fundamental contribution can be 
provided by the data coming from the space telescope 
CoRoT (Convecti on, Rotation and planetary Transits; 
iBaglin et all I 20071 ). This satellite took images of the same 
region of the sky for about five consecutive months reach- 
ing relatively faint stars (up to V~15 mag) in its"exo field", 
with good precision and very high duty cycle. Particularly 
interesting was the first long-term pointing of CoRoT, the 
so-called "Long Run al" (LRal) which imaged a region of 
the sky in the galactic anti-center where a few distant Star 
Forming Regions (SFR) are present. 

With the aim of searching for possible PMS S Set vari- 
ables present in the CoRoT data, we have studied all the 
stars identified a s S Set from the CoRoT variable star clas- 
sification group l|Debosscher et al~1l2009h . looking at near 
infrared (NIR) excess and emission in the Balmer lines, i.e. 
typical featu res of PMS objects. W e mainly used the 2MASS 
photometry l|Skrutskie et al.| [2006) and the yet unpublished 
FLAMESQVLT spectroscopic survey of CoRoT variables 
in LRal performed by C. Neiner and collaborators (ESO 
proposal 082.D-0839A). Among the about 76 stars classi- 
fied as S Set with available spectrum observed by CoRoT 
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Figure 1. Upper panel: light curve of CoRoT 10269979 before 
detrending (see text). Lower panel: selected portion of the light 
curve after detrending with overimposed the least-square fit (solid 
line) to the data obtained with the five "genuine" frequencies (see 

Sect, my 



during LRal, the only candidate PMS 5 Set variable was 
CoRoT 102699796, a relatively faint star with V=15.51 mag, 
located in the direction of Monoceros (RA,DEC 06:43:38.46 
-01:07:46.7, J2000). This paper is devoted to a complete 
characterization of this star on the basis of time-series 
CoRoT observations and a variety of additional data, such as 
VLT mid-resolution spectroscopy and optical-NIR ground- 
based, as well as space infrared photometry. 

In sections 2 and 3 we present the CoRoT photome- 
try and the VLT spectroscopy, respectively. In section 4 
we discuss the evolutionary status of CoRoT 102699796 
and its possible membership to the Star Forming Region 
[FT96] 213. 1-2.2. In section 6 we interpret the observed os- 
cillation frequencies to the light of non-radial pulsation the- 
ory; a summary of the major findings of present work closes 
this paper. 



2 PHOTOMETRIC OBSERVATIONS WITH 
COROT 

CoRoT 102699796 was observed by the CoRoT satellite in 
its "exo field" during the Galactic anti-centre LRal run, 
between October 27, 2007 and March 3, 2008. The data 
span about 131.5 consecutive days. The exposure time was 
of 512 seconds. We used the data corre cted to the N2 
level (the processing steps are described in lAuvergne et all 
I 20091 ). which contains 22036 flux measurements. The first 
step of data pr ocessing consisted in re moving all the points 
with non-0 flag (|Auvergne et al.ll 20091 ). i.e. mainly removing 
those measurements heavily affected by the transit through 
the South Atlantic Anomaly (SAA). A few obvious out- 
liers were rejected, too. The resulting light curve was trans- 
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formed in part-per-thousand (ppt hereafter) using the for- 
mula 1000 x [(Hux/average-1], where average is the mean 
over the whole light curve. This procedure left us with 14997 
useful data points which are shown in Fig. [T] (upper panel). 
The light curve clearly shows a monotonic dimming proba- 
bly due to the ageing of the CCD and/or optics of CoRoT, as 
well as shorter time scale variability, of the order of ~ 15-30 
d and ~3-5 d. The former "periodicity" could be ascribed to 
variable dust obscuration due to the p assing of circumstellar 
mate rial in front of the star (see, e.g. Ivan den Ancker et al.l 
1998), whereas the latter could be due to r otational mod - 
ulation related to coronal activity (see, e.g. ICatalal l2003h . 
However, given that the observed amplitudes are of only 1-2 
hundredths of magnitudes for both types of variability, we 
cannot exclude that all these long term variation are due to 
instrumental effects. Moreover, the light curve also appears 
affected by small "steps" , i.e. sudden jumps in brightness of 
instrumental origin, common a mong the CoRoT li ght curves 
in the "exo field" field (see, e.g. iMislis et alj|2010l ). This oc- 
currence makes it difficult to understand whether a feature 
is intrinsic or due to instrumental problems. To get rid of 
this effect, we decided to apply a high-pass filter to the light 
curve, removing all the long-term variability by filtering out 
all the frequencies lower than 0.50 c/d (all the periods larger 
than 2 days). This very simple approach is effective in re- 
moving also the quoted "steps" because they are not particu- 
larly strong in our case, and no specific treatment is needed. 
We emphasize that the above procedures modify the vari- 
ability features of the star at frequencies much lower than 
those typical of the expected <5 Set pulsation, therefore, they 
do not affect at all the frequency analysis which is presented 
in the next section. 

2.1 Frequency analysis 

The pulsation fre quency analysis was carried out with the 
periodOA package (|Lenz fc Br egcr 2005), which adopts both 
Fourier and least squares algorithms, and permits the si- 
multaneous fitting of multiple sinusoidal variations, and thus 
does not rely on sequential prewhitening. We have also used, 
for comparison purpo ses, the SigSp ec package, which follows 
a different approach (Reegen 2007). 

First, we calculated the spectral window (SW) of the 
data up to the Nyquist frequency (~84 c/d), as shown in 
panel a) of Fig. The S W shows several features due to the 
satellite's orbital period of about 7000 s (variable up to four 
seconds) , and its multiples. This is due to the removal of data 
points during the passage through the SAA. Moreover, since 
the scattered light falling on the detector varies during the 
orbit of the satellite, the periodogram of the star is expected 
to present the orbital period (~13.96 c/d) with the relative 
aliases and harmonics. Frequencies close to these values will 
be rejected in the following analysis. 

The Fourier transform of the data is shown in panel b) 
of Fig. [2] It can be seen that there is a dominant frequency 
at about 6.55 c/d and other signal up to the Nyquist fre- 
quency. Most of the signal is represented by the aliases of 
the dominant frequency and is removed by the prewhitening 
procedure. We extracted frequencies with periodOA using as 
a limit for the last significant frequencies the w idely used 
empirical criteria of S/N=4 (|Breger et al.l Il993h . To esti- 
mate the noise, following well established procedures (see, 



e.g. iRodriguez et all l200rj ; | Garcia Hernandez et alj I2009T ). 
we calculated the average amplitude of the residuals (af- 
ter prewhitening all the significant peaks) in a frequency 
interval of width=5 c/d centred on the corresponding peak. 
Quantitatively the noise resulted to be quite flat beyond 
~1 (c/d) with an average value ~0.03 (ppt). In this way 
we were able to extract 15 significant frequencies. Among 
these, four (at 27.9416, 41.91337 27.9379, and 13.9617 c/d) 
are related to the orbital period of the satellite and will 
be ignored in the analysis. The remaining eleven frequen- 
cies are listed in the first three columns of Table [T] to- 
gether with the amplitude and S/N. The same procedure 
was carried out with SigSpec, obtaining the same results, 
as shown in Table [T] where sig represent the spectral sig- 
nificanc e (see Iftcegen 2007, f or an explana t ion). Following 
iReegenl (|2007f ); IKallinger. Reegen. fc Weissl (|2008h . a value 
of sig—5A6 should be approximately equivalent to S/N=4.0. 
However, as shown in Tab [2 in our case this S/N value cor- 
responds to a sig ~6.6. The origin of such a discrepancy, 
which was already found by other Authors in the literature, 
is not clear. A detailed discussion about this occurrence is 
beyond the scope of present paper, indeed our results are not 
affected at all, since our main discriminant for the significant 
frequencies remains the S/N=4.0 criterion. However, the 
intere sted reader can consult e.g. iGarci'a Hernandez et al.l 
(2009) for a discussion on this subject. 

As a further check on the reliability of the frequency 
extraction, we analysed separately selected portions of the 
time-series, with a size of about 10-15 d, and dealing with 
both "detrended" and raw data (in this case avoiding por- 
tions including "steps"). All the relevant frequencies were 
found in all the cases, reinforcing our confidence about the 
reliability of the results presented here. 

The precision on the single frequencies is also given in 
Table [I] It has been estimated ad opting the defini tion by 
IKallinger. Reegen. fc Weiss] <|2008ft : af = l/(Ty/sig(A), 
where T is the duration of the time-series and sig (A) is the 
significance of the frequency with amplitude A as calculated 
by SigSpec.We also report the uncertainty on the frequen- 
cies on the basis of the Rayleigh criterion (l/4Tj3=0.002 
c/d. This more conservative estimate (at 3-ct level) is the 
value used for the comparison to the theoretical models. 

2.2 Combination of frequencies 

Linear combinations of the terms having the highest am- 
plitudes have been observed in low-amplitud e S Set stars 
both from the ground and from space (see, e.g. iBreger et al.l 
120051 ; iPoretti et al.l 12009) . To search for possible combina- 
tions among the frequencies observed in CoRoT 102699796, 
we decided to use the Combine package, an ad hoc com- 
puter program written by P. Reegen. This program checks 
one frequency after the other for being a linear combina- 
tion of previously examined frequencies. If this attempt fails, 
the corresponding frequency is considered "genuine". Only 
genuine frequenci es are used t o form linear combinations 
subsequently (see Reegen 2010, for a detailed description). 
Applying this package to the investigated star, we find only 

see IKallinger. Reegen. fc Weissl (l2008h and references therein 
for a discussion about this definition 
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Figure 2. Panel a): spectral window of the data up to the Nyquist 
frequency. Panel b): Fourier Transform of the data. Panels c) to 
f): enlargements of Panel b) showing the relevant peaks extracted 
from the periodogram. Labels are the same as in Table [T] 



five genuine frequencies. All the other frequencies can be 
explained as simple linear combination of these terms. We 
varied the input parameters of Combine to verify the ro- 
bustness of the combination calculation. The genuine and 
the combination frequencies are identified in the last col- 
umn of Table [T] and in the periodograms of panels c) to f) 
of Fig. H 



3 SPECTROSCOPIC OBSERVATIONS 

Intermediate-resolution spectroscopy of CoRoT 102699796 
was carried out with .Fib re Large Area Mult i- .Element 
Spectrograph (FLAMES), (|Pasquini et all [2002) mounted 
at the NasmythA platform of the VLT 8.2-m unit tele- 
scope 2. The low-resolution gratings LR02 and LR06 were 
used, leading to a resolving power of R » 6400 and 
R ~ 8600, respectively. The spectral range coverage is 
about 600 A (3960 - 4560 A) for the blue grating and about 
700 A (6440 - 7160 A) for the red grating. 

We acquired four exposures with the LR02 grating, two 
during the night of 2009, February 09 and two others during 
the night 2009, March 22. In order to increase the signal- 
to-noise ratio of the final spectrum, we corrected each sin- 
gle spectrum to the velocity restframe, we verified that no 
residual lines shift is present at our resolution (i. e., due to 
a hidden companion), and finally we combined all of them 
in one spectrum with a total exposure time of 2160 sec, and 
a SNR of ~ 60. 

Only one exposure has been obtained using the LR06 
grating, precisely during the night 2009, February 09, with 
an exposure time of 2400 sec and with a SNR of ~ 40. 

For the following analysis we used the 1-D, wavelength- 



calibrated spectra as reduced by the dedicated Giraffe 
pipelin^E 



3.1 Spectral type and atmospheric parameters 

3.1.1 Determination of effective temperature 

Any attempt devoted to a detailed characterization of the 
chemical abundance pattern in stellar atmospheres relies on 
the accuracy of effective temperature and surface gravity 
determination. 

In this study, we derived the effective temperature by 
using the ionization equilibrium criterion. In practice, we 
adopted as T a g the value that give the same iron abundance 
as computed from a sample of spectral lines (both neutral 
and first ionization stage) present in the spectral range be- 
tween 4000 A and 4500 A (see Tab. 0. 

We used the method of the spectral synthesis in order to 
overpass the strong blending due to relatively high rotational 
velocity ( see below). To perform these calculations we used 
ATLAS9 l|Kurucz 1 1993ft to compute LTE a tmospheric mod- 
els and SYNTHE |Kurucz fc Avrett Il98ll ) to reproduce the 
observed spectrum. Line lists and atomic par ameters used 
in our modeling are fromlKurucz fc Belli lfl995l ) and the sub- 
sequent update bv lCastelli fc Hubrid |2004). From the first 
steps of our iterative procedure, it appeared that the star has 
a low metal content. Thus, we decided to use Opacity Distri- 
bution Functions [ODF] computed for sub-solar metallicity. 
At the end we estimated T cr f = 7000 ±200 K, computed for 
ODF [M/H] = —1. The error on the temperature was esti- 
mated to be the the variation in the parameter that increases 
the x 2 by unity. 

As a by-product, we obtained an estimate of the pro- 
jected rotational velocity usini = 50±5 km/s and an iron 
abundance [Fe/H] = — 1.1 ±0.2. Unfortunately, we were 
able to obtain only an upper limit for the magnesium abun- 
dance [Mg/H] sC -2. 

In Figs. [3] we show two portions of our spectral range: 
one including the Balmer lines and the other covering the 
spectral range 4380 - 4440 A, respectively, with the synthetic 
spectra superimposed. 

It is worth mentioning that the metallicity measured 
for CoRoT 102699796 is rather low with respect to that ex- 
pected on the basis of re cent estimates of the Galactic metal- 
licity gradient (see, e.g. IPedicelli et al.l [20091 . and references 
therein), i.e. [Fe/H]~-0.5. However, as discussed in Sect(5] it 
is possible that the region where CoRoT 102699796 formed, 
is associated with the young open cluster Dolidze 25 and 
with the star forming complex that include other HII re- 
gions (such as Sh 2-284/5), which is known to h ost stars 
as metal poor as [Fe/H]~-0.8 l|Lennon et al.lll990ft . a value 
close to that measured for our target star. The reason of 
such a low metallicity is not clear, and of course needs to be 
investigated in detail, however this is beyond the scope of 
the present paper. 



2 BLDRS vO.5.3, written at the 
http://girbldrs.sourceforge.net 



Geneva Observatory, see 
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Table 1. Pulsation frequencies for CoRoT 102699796. Columns (l)-(4) show frequencies (in c/d and /iHz), amplitudes and S/N derived 
with periodOA, respectively. Columns (5)- (7) prese nt frequencies, amplitudes and siq values obtained with sigSpec, respectively. Column 
(8) show the uncertainty on the frequencies (after iKallinger. Reegen. & ; Weiss (2008)), whereas column(9) exhibit their indentification. 
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Table 3. Physical parameters of CoRoT 102699796 as derived 
from spectroscopy. 
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Figure 3. Top and bottom panels show selected portions of 
our spectrum, i.e., the Balmer lines and the 4380-4440 A 
range, respectively. In both panels the synthetic spectra is over- 
imposed (red line). 



3.1.2 Determination of surface gravity and luminosity 

For early F-type stars, one of the method commonly used 
in the literature for the determination of the luminosity 
class of stars is through the strength of the ratios between 
the spectral blend at AA 4172-4179 A, mostly Fen and Till 
lines, and that at A427 1 A, mostly composed of Fel lines 
|Grav fc Garrison II 19891 ). 

We computed the AA 4172-4179/4271 ratio which re- 
sulted in 1. 60 ±0.15, roughly corre sponding to a luminosity 
class of V (|Grav fe Garrison II 19891 ). 

To derive the surface gravity of our target, we com- 
puted the theoretical behavior of the above ratio as a func- 
tion of log g. After having fixed T e g to the value found in 
Sect. 13.1.11 we computed a grid of ATLAS9 atmospheric 
models with gravities spanning the range between 3.0 and 
4.5 dex. The theoretical equivalent widths and consequently 
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the ratio AA 4172- 4179/4271, have been e valuated using 
the code XLINOP (|Kurucz fc Avrett Ill98ll '). By using this 
curve, we converted our measured ratio and its associated 
error in a measurement of gravity and of its uncertainty, 
obtaining: log g = 3.8 ± 0.4. 

Using the above results, we conclude that 
CoRoT 102699796 has spectral type Fl V. 

The atmospheric parameters derived here can be used 
to estimate the stellar luminosit y. To this aim, we use the 
approach bv lBalona et al.l (| 201 lh . Starting from the Stefan- 
Boltzmann law L=L(T e ff,R) where R is the radius and all 
the quantities are in solar units, we can eliminate the radius 
from the definition of gravity g=g(R,M), where M is the 
mass, obtaining a relation L=L(T e ff ,M,g). Finally, the mass 
can be eliminat e d usin g the empirical mass- luminosity rela- 
tion bv lMalkovl (|2007l) : log L/L =log(M/M ) = 0.00834 + 
0.213 log(L/L Q ) + 0.0107 log(L/L ) 2 , thus deriving the de- 
sired L=L(T e ft,g) relation. 

By using this relation we obtain log(L/Lo) = 1.33±0.3 
for CoRoT 102699796. The relevant stellar parameters for 
this star are shown in Table [3] 

In Figs. [3] we show two portions of our spectral range: 
one including the Balmer lines and the other covering the 
spectral range 4380-4440 A, respectively. The synthetic 
spectra superposed has been calculated using the atmo- 
spheric parameters derived as described before and summa- 
rized in Tab. [3] 
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Table 2. Neutral and first ionization state iron lines identified in the range 4000 - 4500 A for CoRoT 102699796, used for effective 
temperature and iron abundance calculation. For each line we report wavelength and oscillator strength. 
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4 THE EVOLUTIONARY STATUS OF 
COROT 102699796 

In this section we deal with the evolutionary status of our 
target: is it indeed a young intermediate- mass star (Herbig 
Ae type), as we suspect? We have seen in the previous sec- 
tion that the Balmer lines in the spectrum appear to be 
filled. This is the first necessary clue. The next steps are to 
investigate if 1) the star shows infrared excess, and 2) it is 
possible to associate it with a known star forming region. 



4.1 The SED of C0R0T 102699796 

We start by building the Spectral Energy Distribution 
(SED) for C0R0T 102699796. A set of photometric data was 



collected from the literature. Observations in the B, V, R and 
I Joh nson bands are pres ent in the C0R0T EXO-dat cata- 
logue l|Deleuil et al. l2009h . The NIR data in the J, H, and K 
band s were taken from t he 2MASS catalog (Skrutski e et al.l 
120061 : ICutri et al. |[2003 ). Four measurements in the mid- 
and far-infrared were collected from the AKARI Infrared 
Astr onomical Satellite cat alogues IRC (mid-IR all-sky Sur- 
vey, [I shihara et al. 2010T ) and FIS (Far-IR all-sky Survey, 
lYamamura et al. I l2010t ) The data were acquired by the 
Akari-satellite in four bands centered on 9, 18, 65 and 90 /im. 
The position accuracy in the first two bands is estimated to 
be better than 2". At larger wavelengths the pixel size is 
of 26.8" and the average position disagreement between the 
AKARI-FIS catalogue and the Simbad catalogue is of ~6.5 
arcsec. The source identified in the AKARI-FIS catalogue 
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is at a distance of ~23 arcsec from the coordinates given 
in the CoRoT catalogue for CoRoT 102699796. Given the 
pixel sizes of the AKARI camera and the absence of bright 
sources in a radius of 23 arcsec around CoRoT 102699796, 
we can state that the observed FIR flux comes from this 
star. Additional data in the mid-IR can be retrieved from 
the dat abase of the miss ion Midcourse Space Experiment 
(MSX) (|Egan et al. II2003I ). In particular, the target star was 
detected in the A band of the satellite (centered at 8.28 fi), 
which is close in wavelength to the AKARI S9W band. Note 
that these two independent measures are in perfect agree- 
ment one each other (see Fig. 3| . 

The magnitudes of th e star were de r edden ed by us- 
ing the relations given in ICardelli et al. I (|l989l ) between 
the A v and the other photometric bands. In order to es- 
timate A„, we compared the B — V, V — I and V — R 
colours with t h e one s obtained using the stellar models of 
ICastelli et al. I l|2004 ) for the temperature and the metal- 
licity of the star. The E(B-V), E(V-R) and E(V-I) so far 
derived were then converted to A v using the transforma- 
tions of ICardelli et al. I HHI) 

and averaged. The estimated 
extinction is A v = 0.93 ± 0.06 mag. 

The dereddened SED was compared to theoretical mod- 
els of star+disk computed using the software cgplus, a mod- 
eling program, written by CP. Dullem ond, for dusty circum- 
stella r disks , based on the models o f |Chiang fc Goldreich I 
Il997h and |Dullemond et al. I I200II ') . The basic idea is to 
model the disk around Herbig Ae/Be stars with a passive, 
irradiated disk with an inner hole. In our case the model 
which best fit the SED is the one with the parameters re- 
ported in Table [J The SED with the model fit is shown in 
Fig. [4] We note that the luminosity obtained from the in- 
tegration of the star's SED is in excellent agreement with 
that estimated from the spectroscopy, as reported in Tab. [3] 
(see next section for a discussion about the distance of the 
target). 

The shap e of the SED is typ ical of transitional circumstel- 
lar disks (jCalvet et al. I [2005 V These usually have an inner 
region which is depleted of gas and dust and an outer re- 
gion which is optically thick. The strong flux rise around 
10 /J,m is presumably due to the frontal illumination of the 
wall of the opticall y thick outer d isk. F ollowing the classi- 
fication scheme of Sartori et al. I (2010) for Herbig Ae/Be 
stars, CoRoT 102699796 can be identified as a group 2 ob- 
ject, intermediate between embedded objects and stars with 
more evolved dusty disks. These group 2 objects can show 
a doubly peaked SED related to a visible central star (first 
peak) surrounded by a significant amount of circumstellar 
cold matter responsible for the secon d peak in the mid- to 
far-IR band. In the group 2 objects of ISartori et al. I (|2010T l 
the second peak is similar to the first, like in the case of 
CoRoT 102699796 (see Fig.©. 



5 MEMBERSHIP 

CoRoT 102699796 is placed at an angular distance of 1.65 
arcmin from the cente r of the HII region [FT96] 213. 1-2.2 
|Fich fc Terebevlll996l ) which ext ends over 9 arcmin in the 
Monoceros constellation. Moreover I Avedisova I |2002l ) lists 7 
young objects close 4 arcmin) to CoRoT 102699796. We 
believe, indeed, that the CoRoT object, together with the 
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Figure 4. The dereddened SED of CoRoT 102699796 is shown 
with open squares. The continuum line is the total flux given by 
the model star+disk with the parameters presented in Table [4] 



Table 4. Parameters derived by the SED fit. 



Parameter 


value 


Distance 


4000 pc 


Inclination 


60.0° 


Star mass 


2.0 M Q 


Star luminosity 


20.0 Lq 


Inner disk-radius 


4.0 AU 


Outer disk-radius 


400 AU 


Disk mass 


0.06 M 


Inner disk-temperature 


450 K 


Surface density exponent 


-2.0 



7 young objects, belongs to the [FT96]213.1-2.2 HII region, 
which is in turn associated with a big star formation struc- 
ture composed by se veral HII regions, su ch as Sh 2-284, Sh 
2-283 and Sh 2-285 (jRusseil et al. 1120071 ). T he distance es- 
timate of this structure is still controversial. iRusseil et al. I 
|2007l ) derived, using also data from literature, a distance of 
7.89 kpc. Nonetheless one of the biggest part of the star 
formation structure, the HII regi on Sh 2-284, was found 
to be at a distance of 4 kpc by ICusano et al.l (|201ll ). in 
agreemen t within the errors with th e dista nce of 3.6 kpc 
given by iDelgado. Diupvik fc Alfarol <|2010h . Furthermore, 



iRolleston et al. I (|1994| ) also found that two B0V stars in 
Sh 2-285 are at a distance of 4.3 kpc. These results com- 
bined together, place the star formation re gion at a dis- 
tance of 4-4.5 kpc, closer then the estimate of lRusseil et al. I 
|2007h . The luminosity of CoRoT 102699796 derived by the 
frequency analysis (as estimated in Sect 16. 2| ), by the spec- 
troscopic observations and by the SED fitting (see Section 
14. is also compatible with a distance of the object of 4 
kpc. The three colour image shown in Figure [5] obtained by 
combining IRAS images at 25, 60 and 100 /im, illustrates the 
different components of this large star formi ng complex. The 
white open circles are HII regions from the lSharplessl (|l959h 
catalogue, while the green circle-square symbols show young 
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stellar objects from I Avedisova I i|2002l ). CoRoT 102699796 is 
inside the circle of [FT96]213. 1-2.2. 

In the spectroscopic analysis, CoRoT 102699796 was 
found to have a metallicity very similar to that of Sh 2- 
284, supporting the assumption of a common origin of the 
[FT96]213. 1-2.2 and Sh 2-284 HII regions. If confirmed, this 
complex has one of the lowest metallicity galactic known 
SFRs. 



6 INTERPRETATION OF THE OBSERVED 
OSCILLATION FREQUENCIES 

6.1 Calculation of asteroseismic models 

The mo del analysis performed in t his paper is described in 
detail in lDi Criscienzo et al.l (120081) . and was originally sug- 
gested bv lGuenther fc Brownl \2004 ). In summary, given the 
observed N frequencies of the investigated star, we look for 
models having within our grid of oscillation spectra, 

with the following definition of \ 2 



X 



1 \ ^ {fob 



/mod,i) 



(i) 



where f bs,i and fth,i are the observed and the model 
frequency for the ith mode, respectively, characterized by 
an observational cr bs,i an d a theoretical (Tth,i uncertainty. 
In general, the search for the best-fitting model (corre- 
sponding to the oscillation spectrum that minimizes the 
X 2 value) is concentrated in the subgrid of models with 
T e ff and L/Lq within the error box in T e g and L/Lq of 
the star. If the program does not find any model matching 
all the frequencies with x 2 ^li it discards one by one 
the non-matched frequencies, computes the \ 2 f° r the 
remaining frequencies and searches again for models with 
X 2 ^l. Among the resulting models satisfying the condition 
X 2 ^l, the model with the lowest value of the \ 2 ls the 
best-fitting one. The grid of PMS models whe re constructed 
using the ATON code for stellar evolution (|Ventura et al.l 
1 199a) in the standard v ersion for asteroseismic applications 
i D'Antona et alj 120051 ). The physics of the models are 
up-to- date and described in detail in |Pi Criscienzo et al.l 
(2008). The ATON code also allows to incl ude the effects 
of rot ation, according t o the formulat ion bv lEndal fc Sofial 
' 19761) . as described in iMendes et alJ ljl999T l: lLandm et al.l 
2006). The possibility of accounting for the rotation 



effects proves to be extremely useful for this analysis, 
since PMS stars have o n average high rotational velocities 
(Bo hm fc Catalall 19951) . In general, what we may determine 
is the projected velocity vsini, but in this particular case we 
can derive the true rotational velocity (~60km/s) from the 
inclination of the disk, assuming the disk is in the stellar 
rotation plane, and this can be used to further fix the best 

fit model. 

We stress here that the approach bv lEndal fc Sofial (|l976T l 
accounts only for the hydrostatic effects of rotation, ne- 
glecting the internal angular momentum redistribution; 
in this work we use a rigid body rotation. The initial 
angular momentum of the star is provided as a physical 
input, and chosen to reproduce the observed surface angu- 
lar velocity of the star at a given position in the HR diagram. 
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Figure 6. HR diagram showing the spectroscopic position of 
CoRoT 102699796 (open circle) with the respect to several PMS 
rotational tracks taken from the model grid. The black full cir- 
cles mark the position of the best fit models found comparing 
observed and predicted frequencies with m=0 while magenta tri- 
angles are the models described in Table 6. The red solid lines 
and the green dashed lines show the 5 Set and 7 Dor instabilit y 
strips by iBreger &: Pamvatnvkhl jl998h and iGuzik et al. 
respectively. 



Models were computed from 1.5 to 2.3 Mq, because 
these are the masses that enter in the empirical box in the 
HR diagram (see Fig.©, with a step of AM=0.01 M . The 
evolution begins on the Hayashi track, and the calculation 
ends when the star has reached the MS. We note that ne- 
glecting the protostellar evolution and its effect on the inner 
structure and pulsation properties does not change the con- 
clusions of the paper because the investigated star lie on 
the horizontal portion of the PMS evolutionary tracks, sub- 
sequent to the thermal relaxing phase undergone by stars 
around 2 Mq in the earliest phases of the PMS evolution. 
The step in effective temperature between consecutive mod- 
els of the grid is AT cS ~50K. 

We use the LOSC oscillation code bv lScuflaire et all (|2008l ) 
to calculate the frequencies of radial and non radial oscilla- 
tion modes for each evolutionary structure of the grid. We 
limit the spherical degree to 1<4 and contrary to previous 
works on PMS pulsating stars we extend our pulsational 
analysis also to gravity modes since the oscillation spectra 
of this star has an anomalous extension to low frequencies 
suggesting that at least the lowest ones could be low order 
g-modes. 
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Figu re 5 . Com position of three IRAS image at 25, 60 and 100 fim. White open circles show the HII region s acco rding to the 
ISharpless I Jl959h catalogue, while green circle-square symbols show young stellar objects reported by lAvedisova I (120021 ') . The target 



(separations lower than 0.7/xHz or 0.06 c/d) present in 
the power spectra of the majority of the we ll-studied 5 
Set stars, as noted by iBreger fe Bischof I (|2002l ) . Since the 
actual cause for this phenomenon is still unknown we prefer 
to remove from our fitting procedure. 
In this way we have four observed frequencies to be fitted 
with computed ones in order to find the best fit model 
through the \ 2 analysis. 



We assume as uncertainties on observed frequencies in 
eq. [T] values which are 3-ct the observational errors given in 
Section 2.1 and consider conservatively theoretical errors 
of the same order of magnitude of the observational ones, 
taking into acco unt all the in trinsic numerical uncertainties 
of models (see iMoval 120071 . for a detailed quantitative 
description) . 

As a result we find that the best fit model corresponds 
to a very young star (t~3 Myr, with no convective core 
yet) with M=1.74 M Q , T to ff=6850K (CASE A in Table ©. 
We note that the % 2 reported in table is calculated on the 
basis of the first three frequencies because unfortunately this 
model does not match frequency f5, which however is the 
one with the lowest S/N ratio (see Tab.[TJ. Moreover, if we 
allow x 2 to be > 1, we find that a much more massive model 
(described in Table [5] as CASE B) minimizes \ 2 ■ According 
to this model, only two frequencies are pure g-mode of low 
order. 

We remark that in the analysis described above the ro- 
tation of the star (see Sect. 3) was considered only in the 
computation of the evolutionary model. For completeness we 
recall that if / denotes the frequency in the absence of ro- 
tation, a slow solid rotation with angular velocity Q slightly 
alters the frequency with 1>0 in the following way: 



star CoRoT 102699796 is inside the circle of [FT96]213.1-2.2. 
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Figure 7. Comparison between the observed spectrum of star 
CoRoT 102699796 and the spectrum of the best fit model ob- 
tained in the two cases indicated in Table[5] Filled circles and tri- 
angles, crosses and asterisks indicate theoretical frequencies with 
1=0, 1, 2, 3, respectively. The open circles indicate the matched 
frequencies 



6.2 Comparing observed and predicted 
frequencies 

Among the observed frequencies reported in Table [T] we 
try to fit those which are not consistent with a linear 
combination of the others, namely fl to f5, in the order 
adopted in Table [T] We also start excluding from the 
analysis fa which is too similar to /i and we suggest that 
these two frequencies are an example of frequency pairs 
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fl,n,m = fi,n + m/3i yn Q (2) 

where (3 is a qu antity described in detail in 
|Pi Criscienzo et al,l (|2008h and calculated by the adiabatic 
oscillation code for each pair (l,n) . We then relaxed the pre- 
vious constraint used to find the best fit model by adding 
the multiplet (triplet and quintuplet, etc. according to eq.[2]) 
frequencies to our predictions, and we compared again the- 
ory and observations using the same technique as above (but 
this time trying to fit also fs). As a result, the number of 
best-fitting models grows in number, as shown in Tab. [(J 
where we list all the models between Xmin an d xLn + 1 (i- e - 
within 1 a). However, when we report these models in the 
HR diagram (empty triangles in Fig. [6]), we find out the very 
interesting result that these models cluster around or close 
to the two best fit models previously found, namely CASE 
A and B. This occurrence suggests that models consistently 
predict two possible solutions for the position of the star in 
the HR diagram: a) with LogL/L Q ~ 1.24-1.3 and M/M ~ 
1.8; b) with LogL/L© ~ 1.6 and M/M Q ~ 2.3. To decide 
which "class" of models has to be preferred, we note that 
the "class a)" are the favoured models because: i) they pro- 
duce the lowest values for the \ 2 an d ii) their luminosities 
are in better agreement with the values inferred from spec- 
troscopy and from the integration of the SED (although at 
least formally also the luminosity of class b) models are in 
agreement with the spectroscopic value at 1 a). On the ba- 
sis of the above considerations, our best-fit model is the one 
that minimize the x 2 (model al in Tab. [6]), with mass, ef- 
fective temperature and luminosity of 1.84 M/Mq, 6900 K, 
and 1.29 LogL/L©, respectively. Similarly to CASE A, it is 
a young model (~2.5 Myr) still without a convective core. It 
is important to note that the observed frequency spectrum 
can be interpreted in terms of g-modes with low-moderate 
n- value for class a) models, whereas class b) models are char- 
acterized by both p- and (/-modes. 

These results, if confirmed, would be of crucial im- 
portance because it would be the first time that gr-modes 
are detected in a Herbig Ae pulsating star. An even 
more complicate picture emerges if we consider that 
CoRoT 102699796 lies (see Fig. [6]) in the region of the HR 
diagram where S Set and 7 Dor inst ability strips intersect 
(at le ast for solar metallicity, see e.g. lHandler fe Shobbrookl 
2002). We recall that the frequency spectrum of 7 Dor stars 
is interpreted in terms of high radial order (/-modes. Our 
findings suggest that C0R0T 102699796 shows a frequency 
spectrum with intermediate properties between S Set and 
7 Dor, characterized by frequencies lower and higher than 
those typical of 5 Set and 7 Dor variables, respectively. In 
fact, according to our interpretation this star is a pre-main 
sequence pulsator which shows gr-modes of low-moderate 
radial number and order in its frequency spectrum (or both 
p- and (/-modes, see above). In this context, it is worth 
noticing that hybrid S Sct-7 Dor pulsators were already 
found among more evolved stars by seve r al authors(see, e.g . 
Henry fc Fekel 120051 : iKing et al l l200fj : iRowe et all 120061 : 
Handler! |200S| ). Moreover, the recent data obtained with 
the KEPLER satellite seem to show tha t the hybrid phe- 
nomenon is more common t han expected (|Grigahcene et all 
l20ld : ICatanzaro et al.ll2010l ). 

However, prior to arrive to firm conclusions concerning the 
nature of C0R0T 102699796, a more detailed theoretical 



study is necessary by using a non-adiabatic oscillation 
code for a stability analysis of the selected modes. Fur- 
thermore the main f eature of non-adiabatic codes as MAD 
l|Dupret et al. Il2003l ) is their ability to determine accurately 
the amplitudes and phases of variation of different physical 
quantities in the very outer layers of the stars. In particular 
the normalized amplitude and phase of effective tempera- 
ture variations can be determined, significantly improving 
the accuracy of the mode identification methods through 
multi-band photometric approaches. An important point is 
that the non-adiabatic predictions are very sensitive to the 
modelling of the thin convective envelope of these stars and 
to the treatment of its interaction with pulsation, thus only 
time dependent convective versions of these codes would be 
able to reproduce both the blue and the red edge of the 
instability region of each mode. 



7 DISCUSSION AND CONCLUSIONS 

In this paper we presented a comprehensive study of 
C0R0T 102699796, a star observed by the C0R0T satellite 
during its first Long Run in the anticenter direction. The 
star was classified as 5 Set and shows anomalous Balmer 
lines at a first analysis of mid-resolution spectra obtained 
with GIRAFFE@VLT, thus being a likely candidate of the 
PMS S Set class. 

We analysed the time-series photometry observed by 
the satellite C0R0T, detecting the presence of five indepen- 
dent oscillation frequencies in the range 3-6.5 c/d. These low 
values are somewhat atypical for PMS and "normal" 5 Set 
variables, whose peaks are usually placed around 10-30 c/d. 

The star's physical parameters were estimated on the 
basis of the quoted GIRAFFE® VLT mid-resolution spectra, 
finding that the star has spectral type F1V, T cff =7000±200 
K, log(g)=3.8 ± 0.4, [M/H]=-l.l±0.2, usini=50 ± 5 km/s, 
L/L =21±^. Therefore C0R0T 102699796 is the first 
intermediate- mass PMS pulsating star for which poor metal- 
licity has been reported. 

The evolutionary status of the target star has been ex- 
plored in detail by constructing the SED, extending from 
the optical to the mid-infrared, using ground-based and 
satellite data. As a result, we find that the shape of the 
SED is typical of transitiona l circumstellar disks. Following 
the classification scheme of ISartori et al. I (|2010T ) for Her- 
big Ae/Be stars, C0R0T 102699796 can be identified as a 
group 2 object, i.e. between embedded objects and stars with 
more evolved dusty disks. Its PMS nature is reinforced by 
the possible association with the HII region [FT96]213. 1-2.2 
(|Fich fc Terebevlll996l ). 

The pulsation frequencies have been interpreted in the 
light of t he non-radial pulsati on theory, using the LOSC 
program (|Scuflaire et al.l l2008h in conjunction with static 
and rotational evolutionary tracks cal culated by means o f 
the ATON code for stellar evolution (| Ventura et all 1 19981 ) 
in the standard version for asteroseism ic applications with 
low metallicity l|D'Antona et al.ll2005l ). A minimization al- 
gorithm was used to find the model that matches the ob- 
served frequencies. As a result, either including or not the 
rotation, we find two possible solutions for the position of 
the star in the HR diagram: a) with LogL/L© ~ 1.24-1.3 
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Table 5. Mode identification considering all frequencies with m=0 













f obs,l 


f obs,2 


fobs, 4 


f obs,5 












75.91 


38.40 


78.62 


41.46 












(MHz) 


(/.Hz) 




(MHz) 


CASE 


x 2 


M/M 


T cff (K) 


LogL/L 


/tfe,l(l,n) 


/th,a(l.n) 


/«i,4(l,n) 


/th,s(l,n) 


A 


0.83 


1.74 


6850 


1.21 


75.89(1,-2) 


38.36(2,-9) 


78.67(3,-4) 


41.98^(1,-4) 


B 


4.92 


2.29 


6950 


1.59 


75.66(3,0) 


38.57(3,-6) 


78.59(2,0) 


41.33(2,-3) 



Table 6. Mode identification adding in the simulated spectra multiplet frequencies 













fobs,l 


fobs,2 


fobs, 3 


fobs, 4 


fobs, 5 












75.91 


38.40 


75.61 


78.62 


41.46 












(/.Hz) 


(MHz) 


(MHz) 


(MHz) 


(/.Hz) 


MODEL 


x 2 


M/M 


T cff (K) 


LogL/L 


l,m,n 


l,m,n 


l,m,n 


l,m,n 


l,m,n 


al 


0.28 


1.84 


6900 


1.29 


3, 1,-3 


3,1,-8 


3,2,-5 


3,-2,-2 


3,-1,-8 


bl 


0.37 


2.28 


7050 


1.59 


3,-1,0 


3,1,-7 


2,-2,0 


2,-2,0 


2,2,-5 


b2 


0.46 


2.36 


7000 


1.63 


3,-2,1 


3,2,-7 


2,0,0 


2,1,0 


2,2,3 


a2 


0.84 


1.90 


6950 


1.34 


3,3,-5 


1,0,-3 


2,1,-2 


3,-1,-2 


3,2,-4 


b3 


0.86 


2.32 


7250 


1.62 


2,1,0 


2,0,-3 


3,2,0 


2,2,0 


3,-1,-3 


a3 


0.88 


1.83 


6875 


1.28 


3,1,-3 


1,-1,-3 


2,0,-2 


3,-2,-2 


3,2,-15 


bl 


0.90 


2.35 


6850 


1.61 


2,1,0 


2,0,-3 


3,2,0 


2,2,0 


3,2,-6 


b5 


0.90 


2.38 


6950 


1.64 


3,-1,1 


2,0,-3 


3,2,0 


2,2,0 


2,2,-6 


b6 


1.00 


2.36 


6900 


1.61 


3,-1,1 


2,0,-3 


3,2,0 


2,2,0 


2,2,-6 


al 


1.14 


1.81 


6975 


1.27 


2,-1,-2 


1,-1,-4 


3,-3,-5 


3,-3,- 2 


3,3,-19 



and M/M ~ 1.8; b) with LogL/L ~ 1.6 and M/M ~ 
2.3. Models with lower mass and luminosity are preferred 
because: i) they minimize the \ 2 an d ii) their luminosities 
are in better agreement with the values inferred from the 
spectroscopically estimated logg and from the integration of 
the star's SED. Finally, the model named al in Tab.[B]is our 
best-fit model (minimum % 2 using all the five significant fre- 
quencies), having mass, effective temperature, and luminos- 
ity equal to 1.84 M/M , 6900 K, and 1.29 LogL/L , respec- 
tively. The position in the HR diagram of this model is in 
excellent agreement with the spectroscopic one. The model 
is very young (age~2.5Myr) without a convective core, in 
agreement with the shape of the SED that shows a signif- 
icant FIR excess. The theory- observation match is based 
on the interpretation of frequencies /i -5- fs as g-modes of 
low-moderate n-value. To our knowledge, this is a new re- 
sult and the first time that such modes are observed in PMS 
intermediate-mass pulsating stars. However, this is not com- 
pletely surprising since it is also the first time that a so low 
frequency range is observed in a pulsating Herbig Ae star. 
In this context, it is also remarkable that CoRoT 102699796 
lies in the region of HR diagram where 8 Set and 7 Dor 
instability strips intersect. Hence we are probably dealing 
with a star showing intermediate (hybrid) characteristics in 
between the two variability classes. However, we note that 
if the solution predicted by class a) models applied, the star 
could not properly be defined as an hybrid because only 
g-modes would be present. In this case it would be better 
defined as a PMS anomalous 7 Dor star. 
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